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Ugitech S.A 



UGITECH, a main industrial stainless steel long 
products actor based on a famous legacy 



Merger of Usinor, 
Arbed and Aceralia 
to create Arcelor 



Schmolz & Bickenbach 
acquires Ugitech 

i 

2006 I 



Foundation of Sprint Metal 
through the merger of the wire 
drawing Facilities of Ugine 
Aciers and Creusot-Loire 



Ugine-Savoie 
acquires Bedini 

i 



First heat with the 
Ugine Perrin Process 

1 

1925 I 





Ugitech creation to merger 
of Ugine-Savoie Imphy and 
Sprint Metal 



New dust cleaning 
Machine in Ugine 



1972 

T 



Merger of Usinor and Sacilor 
Ugine Aciers is renamed 
Ugine-Savoie 



1908 

T 

Paul Girod begins the 
production of alloyed 
steels at Ugine 



Merger of Ugine the 
Pechiney Group to form 
Ugine Aciers 
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Dedicated to stainless steels long products 

UGITECH Operations 



f 



Production / R&D 



UGITECH (Ugine) 

• Melting Shop 

• Remelting Furnace 




Transformation 



Bar Drawing 

• UGITECH (Ugine) 

• Trafilerie Bedini (Milan) 



II 



Wire Drawinq 




• UGITECH (Bourg-en- 


Bresse) | 


• UGITECH (Brionne) 


1 1 


• UGITECH (Imphy) 


M 


• SMEZ (Reichshof) 




• SMEZ (Brumby) 





Distribution + Services 



SCHMOLZ + BICKENBACH Network 



UGITECH (Grigny - Rouen - Marignier)| 
UGITECH Sri (Milan) | | 

UGITECH GmbH (Renningen) 
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A Heat Transfer Model for Ugitech's Continuous Casting Machine 



Ugitech: 

- Member of the S+B Group 



CHAUFTAGE TRANSLATION 
REPARTITELJR R£PART|TEUR 

— Dedicated to Stainless Steels Long Products 



OSCULATION IINGOTIERE 



The Continuous Casting Machine: 

- Technological features 

- Problems to be solved 

Til ■ 

Numerical Simulation: 

- The heat transfer model tG600 

- Some results / comparisons with measurements 

- Recent developments : mechanical model 



TETE MAN 



6 200 



Conclusion: 

- Work in constant progress 

- Generalization of numerical simulation 



ENSEMBLE MA 



4 * 




STOCKAGE 
Fr DISTRIBUTION 
AUTOMAT1QUE 
DE POUDRf 
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The Continuous Casting Machine: 



Liquid metal 



Technological features: 

- A 3-strand fully vertical machine 

- Bloom 205 2 mm with 040 mm rounded corners 

- Casting of stainless steels (any grade) 

- Casting speed around 1 m/min 

- Lubrication with mineral mould flux 

- Metallurgical length: 14 m 

Phenomena 

- Casting of superheated liquid steel: +40°C (=1400°C) 

- First solidification within the mould, through a thin layer of 
infiltrated molten mould flux 

- Passing through several water cooling nozzles 

- End of solidification through radiation mainly 

The blooms are torch-cut after 14 m, in lengths about 3,4 m. 
This constitutes a maximum limit we cannot exceed as far 
as the end of solidification is concerned. 



Water cooled mould 
length 800 mm 



Flat nozzle cooling 



Round nozzle 
cooling 



Holding rolls 



Round nozzle 
cooling 



Natural cooling 



Liquid well 
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Problems to be solved i: 



Physical phenomena 

- Turbulent fluid flow before solidification 

- Heat transfer, with or without convection 

- Mechanical stresses and strains after solidification 

As a matter of fact, all these phenomena are tightly coupled and 
extremely non-linear. Solidification and mechanical contacts 
(mould/rollers) are a modelling difficulty. 

Model simplifications 

- The fluid flow is calculated separately to check the vertical extend 
of the radial convection. Only a local increase in the heat 
conductivity brings the convection influence in the model 

• Only the upper part of the mould content is affected by convection 

- All the heat exchanges are in the horizontal plane 

• This is not true in the upper part... But we verify that the overall heat 
extracted through the mould is correct, compared with measurements 

- All the stresses and strains are in the horizontal plane 

• Fairly correct if one considers that things change slowly along the z 
axis. 
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A Heat Transfer Model for Ugitech's Continuous Casting Machine 



Ugitech: 

- Member of the S+B Group 

- Dedicated to Stainless Steels Long Products 



The Continuous Casting Machine: 

- Technological features 

- Problems to be solved 

Numerical Simulation: 

- The heat transfer model 

- Some results / comparisons with measurements 

- Recent developments : mechanical model 

Conclusion: 

- Work in constant progress 

- Generalization of numerical simulation 



mould 




ring 



SCI 



SC2 



l | # 
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The heat transfer model (1) 



I20mri 



Symmetry 



205mm 



Symmetry 



Symmetry 
Symmetry 




A 2D model : geometry and meshing 

- Considering the geometry of the problem, the heat is supposed to 
flow within the section of the bloom 

■ We neglect the 3D aspects. . . 

• True along the length of the product 

• To be corrected near the top meniscus 

■ The symmetries are used whenever possible 

• Only one eighth of the section is meshed 

• Whenever required, we can use a full section mesh 




Metal domain, with a fine regular 
mesh, to keep track of the 
solidification front from 100% 
liquid to lOOu. solid 
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The heat transfer model (2) 



Domain settings 

- The solidifying stainless steels undergo different phase changes, often involving more than one 
solid phase (5 and a ferrite, y austenite) 

- The solidification interval is modelled using equilibrium thermo-dynamic calculation. A solid 
fraction is used to calculate the thermo-physical data as a composite between liquid, and solid. 

- For the time being, the local solidification time is not taken into account 



... but it's part of our projects... Or thermo-chemistry software allows for it. 



T°C solidus 



T_liq 

T_sol 

DelT 

Tadim 

cp_Sc 

cp_Lc 



The curves are fed in Comsol Multiphysics either as: 

• Direct scalar expressions 

• Formal derivative against T 

• Tables 

• or a mixture of all 



1456 . 06 [degC] 
1407 . 76 [degC] 
T_liq-T_sol 



(T-T_sol) /DelT 

3197 . 685121 [J/ (kg*K) ] 

107 . 528668 [J/ (kg*K) ] 



Adimensional temperature 
Solid heat capacity (order 0) 
Liquid heat capacity (order 0) 



Liquidus temperature 
Solidus temperature 
Solidification range 




T°C liquidus 



Adimensional solid fraction as a function 
of the adimensional temperature 
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The heat transfer model (3) 



Some thermo-physical data 




H s =C Psc (T-T„,) + Cp 



SP 



rp% rr< 2 rri 3 r T~' 3 

1 ~ ' sol , 1 ~ 1 sol 



Heat conductivity W/m/K 



calc 

— (k SP T + k sc )f fra + k 

mouvL 

((k LP T + k LC )(l-f fm )) 



k solid 




Solid fraction 




k liquid 



kmouvL = 1 + (MouvL - \)flc2hs{M^ m - Position, M liss ) 
Correction factor for convection effects 



2 



Sp2 



H L = Cp LC (T - T lia ) + Cp LP — ^ + dHsol 



H calc =f fra H s+ (l-f fra )H L 
Cp(T) 



dT 



Heat capacity J/kg/K 




1400 
Temperature °C 



November 6 th 2008 



Comsol Conference Hannover 2008 



11 



Providing swciei steel bo jtiors 



The heat transfer model (4) 



Boundary conditions 

- Depending on time 

• Since the simulation is 2D, the bloom section is translated within the 
machine, at the casting speed... The boundary conditions vary in time, 
to rebuilt the history of the section. 

- Four kinds 

• Conduction in contact with the upper copper mould 

- A simple extracted flux equation depending on the z/t axis 

• Radiation along the free length 

- A temperature dependant equation 



Conduction 



Convection 




0.95 

0.9 ^ 
0.85 

0.8 - 
0.75 - 

0.7 
0.65 ^ 

0.6 

0 200 400 600 800 1000 1200 1400 

Surface emissivity as a function of surface temperature 

Volumic heat flux 

- Flux generated by an electromagnetic stirring device 




EMS 



Radiation 
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UGITECH 
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• Convection: 

- The capacity of Comsol to handle large equations allows us to describe each cooling 
nozzle independently. The (huge) equation used, gives the value of the heat transfer 
coefficient as a function of x (horizontal coordinate on the bloom surface), z/t (vertical 
position), the water flow rate, and the geometrical characteristics of the nozzle. Each 
nozzle has its own equation... 



hLocal=157 0* (((((((( (x+ (Nspray-1) /2*dSpray) / (Dist . *tan (AngH*pi/360) ) ) . A 2 + (t/ (Dist . *tan (AngV*pi/360 ) ) ) A 2) <=1) ) . * (cos (at an 
(sqrt ( ( (x- (- (Nspray-1) 12) *dSpray) . A 2) +t A 2) /Dist) ) ) + (Nspray==2) *((((( (x- (1- (Nspray-1) 12) *dSpray) / (Dist . *tan (AngH*pi/3 60 ) ) 
) . A 2+ (t/ (Dist .*tan (AngV*pi/3 60 ) ) ) A 2) <=1) ) . * (cos (at an (sqrt ( ( (x- ( (Nspray-1) 12) *dSpray) . A 2) +t A 2) /Dist) ) ) ) ) / (2*pi* (Dist A 2) *N 
spray* (Dist . *tan (AngV*pi/360) ) / (Dist . *tan (AngH*pi/360) ) * (sin (AngH*pi/3 60 ) *tan (AngH*pi/360) + (cos (AngH*pi/3 60 ) -1) ) ) *Qeau*N 
spray* 0 . 06) )/3.6) . A 0.55*(l-0.0075* (Teau-2 73 . 15) ) /4 




Round nozzle Paired round nozzles Flat nozzle 



Values of h - convection coefficient -- as a function of x and z 
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Some results ... (1) 



iime=2!9. 997531 Surface: Temperature Isovaleur: Temperature 

J ■! 



M.rai: !.!OT 
1500 




Maxi: 1475 
1475 



Solidification 
front 



Mini: 950 IlKfSffi 



Time=219. 997531 Surface: Flux de chaleur total [Wj'nf 
iicivaleur: Plus de chaleur total [W/rn 2 ] 



0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 
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Mould 




Ring 



Mff1lf**Me-3 Mini:2.50e4 



Temperature profile in an 8 th of a bloom section 



Heat flux profile in an 8 th of a bloom section 



Roller mark 




Grade 304L 
Casting speed 1,12m/min 



Section displayed > | 



The calculation time is about 1000s 



Natural cooling 
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Some results ... (3) 



• The length of the liquid well (end of solidification at the axis) can be calculated: any 
torch cutting on the liquid is avoided. 



• The surface defects are also correlated with the temperature gradient and evolution 




1 ■ 1 ■ • ■ ] J ^ Metallurgical length [m] 

0 2 4 6 8 10 12 14 16 

Position [m] 



Temperatures along the strand at the corner, at mid-face and at this axis 
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UGITECH 
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... and comparisons with measurements (1) 

• Pyrometric measures have been taken along the CC strand and fall within ± 10°C from the forecast. The surface 
temperature is well calculated 

1500 t- —I 




900 H 1 1 1 1 1 1 1 1 1 — 

0 2 4 6 8 10 12 14 16 18 

Metallurgical length (m) 
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... and comparisons with measurements (2) 



An incidental torch cutting displaying a liquid metal drop, confronted with the casting conditions shows that the liquid 
well position is also correct. 



1500 
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Recent developments : thermo-mechanical model (1) 

The most difficult part of the model resides in the 
mould, where one creates the solid phase, and where 
all the phenomena occur: 

- Fluid flow (which changes the heat fluxes locally) 

- Phase change which induces volume changes 

- Variable, non-linear surface contact between mould and 
metal 





Most of the defects are coming from this early 
solidification 

- Cracks 

- Segregations near the product's skin 

- Depressions 

- Oscillations marks 

If not to give an exact forecast, at least does a model 
assist the researcher in understanding the process 





J) 

I - 
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Thermo-mechanical model (2) 



Two dialoging meshes: 



Convection (water) 



Contact condition 



Rair — 6air/A,air 



Rl — 6s/A-laitier.solide ^2 — S|/^-laitier.liquide 

Rair = (T-£laitier-£metal'(T|aitier+3.T| a jtier-Tsurf+3.T|aitier-Tsurf+Tsurf) 




Ferrostatic pressure 
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• Material : 

- Elasto-plastic material, with a perfect plastic behaviour in the solid phase 

• This is a theme for our current work... A visco-plastic behaviour would be better... 

- Purely elastic behaviour, with a very low Young modulus in the liquid phase 

• This does not modify the results since the Yong modulus chosen is several orders of 
magnitude lower than in the solid phase 

• This assumption avoids strain accumulation in the liquid, while transmitting correctly the 
pressure at the solidification front 

• Boundary conditions 

- Symmetries (conventional BC) 

- Convection coefficient at the external boundary of the mould. 

- Ferrostatic pressure at the inner boundary of the mechanical mesh 

- Thermal resistance between mould and metal 

• Taking into account the mould flux properties and thickness 

• Taking into account the eventual air gap deduced from the deformations of the surface 

- Contact condition at the surface of the metal 

• Hand built to accommodate the transient calculation 

• Based on a smoothed reaction force proportional to the ferrostatic pressure 
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Thermo-mechanical model (4) 
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Air Gap 



Heat Flux 



Temperature 



Meniscus 



First solid 



Mould 




Air gaps 





1 .4mm 



400kW/m 2 



1220^ 



November 6 th 2008 



Comsol Conference Hannover 2008 



22 



Prouiainf) s»ci6l sleet BOljTkxis 



Thermo-mechanical model (5) 



Results: 



Flux de ehaleut extrait ■->:■ fan. 




Heat fluxes along the length of the mould 

These results are qualitative for the time 
being. The big difficulty of convergence and 
numerical stability has been overcome and 
work is still in progress regarding the 
constitutive laws of the near melting solid 

The calculation time is about 5 hours 



Temperature and 
iso-solid fractions 



Stresses and strains 




Contact zone 



Solidification front 



Iso solid fraction 




Air gap 



Section of the product at the mould exit 
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A Heat Transfer Model for Ugitech's Continuous Casting Machine 

Ugitech: 

- Member of the S+B Group 

- Dedicated to Stainless Steels Long Products 

The Continuous Casting Machine: 

- Technological features 

- Problems to be solved 

Numerical Simulation: 

- The heat transfer model 

- Some results / comparisons with measurements 

- Recent developments : mechanical model 

Conclusion: 

- Work in constant progress 

- Generalization of numerical simulation 
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• Using our current thermodynamical knowledge, and Comsol: 

- We can calculate the heat transfer of any heat's solidification 

- Our heat transfer model is accurate enough to track the liquid well 

• It has been checked against real data on two grade families 

- We have first results regarding the interfacial heat transfers in the mould 

• Constant progress on three battle fronts: 

- Apply the use of the model to more and more case studies 

• Casting speed limits, 

• Format changes 

• Asymmetry 

• New grade families 

- Feed the model(s) with accurate data (thermo-physical data acquisition) 

• Increase our current databases further than 304, 326, 420, etc. 

• Measure p, Cp, X, solidification enthalpy for new grades 

• Improve the nozzles description 

- Improve the model itself 

• constitutive laws for near melting steels 

• 3D ?? 
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Conclusion (2) 



Considering the success of this model, we are extending our numerical 
simulation effort to other fields: 



Surface: Ttot \°C] 



- Ladle simulation (2D-axi heat transfer, 3D-CFD) 

- Localised stresses on a machined part (3D elastic) 

- ESR heating (3D AC-DC + heat transfer) 

- Bloom cooling after CC (3D heat transfer) 

- Tundish simulation (3D k-e CFD) 




November 6 th 2008 



Comsol Conference Hannover 2008 



26 



